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Abstract 

We will discuss the spectrum of the eta mesons making use of the Nambu-Jona-Lasinio (NJL) 
model supplemented with a model of confinement. We will go on to discuss the properties of 
mesons at finite temperature and the phenomenon of deconfinement. We will then discuss some 
excited states of the quark-gluon plasma calculated in lattice QCD models. These resonances are 
thought to be created in heavy-ion collisions. We consider the role these states play in leading to a 
hydrodynamic description of the plasma at early stages of its formation. 
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I. INTRODUCTION 



In this presentation we will discuss the properties of QCD as one moves upward along 
the temperature axis from the point at T = 0, where the matter density is zero. Note that 
for experiments at RHIC the associated chemical potential is small. 
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Fig. 1 Schematic phase diagram of nuclear matter (Taken from F. Karsch, Lect. Notes 
Phys. 583, 209 (2002) [hep-lat/01060T9| .) 

II. PROPERTIES OF THE 77 MESONS CALCULATED WITH THE NJL INTER- 
ACTION AND A MODEL OF CONFINEMENT 



C. M. Shakin and Huangsheng Wang, Physical Review D, 65, 094003 (2002) 



Recent work has shown that the singlet-octet mixing angles of the 7/(547) and r/(958) 
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are different. That may be demonstrated either in extended chiral perturbation theory 
or by analysis of a large body of experimental data. The conclusion is that the ?7(547) 
is almost entirely of octet character, while the 7/(958) is mainly of singlet character with 
about 10% octet component. It is possible to calculate the mixing angles and decay constants 
in our generalized Nambu-Jona-Lasinio (NJL) model, which includes a covariant model of 
confinement. Our model is able to give a good account of the mass values of the 77(547), 
t/(958), 77(1295), and 7/(1440) mesons. (We also provide predictions for the mass values 
of a large number of radially excited states.) It is well known that the Ua(X) symmetry 
is broken, so that we only have eight pseudo Goldstone bosons, rather than the nine we 
would have otherwise. In the NJL model that feature may be introduced by including 
the 't Hooft interaction in the Lagrangian. That interaction reduces the energy of the 
octet state somewhat and significantly increases the energy of the singlet state, making it 
possible to fit the mass values of the 7/(547) and t/(958) in the NJL model when 
the 't Hooft interaction is included. In this work, we derive the equations of a covariant 
random phase approximation that may be used to study the nonet of pseudoscalar mesons. 
We demonstrate that a consistent treatment of the 't Hooft interaction leads to excellent 
results for the singlet-octet mixing angles. (The values obtained for the singlet and octet 
decay constants are also quite satisfactory.) It may be seen that the difference between the up 
(or down) constituent quark mass and the strange quark mass induces singlet-octet mixing 
that is too large. However, the 't Hooft interaction contains singlet-octet coupling that 
enters into the theory with a sign opposite to that of the term arising from the difference of 
the quark mass values, leading to quite satisfactory results. In this work we present the wave 
function amplitudes for a number of states of the eta mesons. (The inclusion of pseudoscalar 
axial-vector coupling is important for our analysis and results in the need to specify eight 
wave function amplitudes for each state of the eta mesons.) We present the values of the 
various constants that parameterize our generalized NJL model and which give satisfactory 
values of the eta meson masses, decay constants, and mixing angles. It is found that the 
calculated mass values for the 7/(1295) and 7/(1440) are quite insensitive to variation of the 
parameters of the model whose values have largely been fixed in our earlier studies of other 
light mesons. 
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III. CHIRAL SYMMETRY RESTORATION AND DECONFINEMENT OF 
LIGHT MESONS AT FINITE TEMPERATURE 



Hu Li and C. M. Shakin , \hep-ph/0209136 



Confinement model: 



V (r) = nr exp(— fj, r) , 



(3.1) 



V°(k — k') = -8tvk 



1 



\(k -k') 2 + /i 2 ] 2 [(k — k') 2 + /i 2 ] 3 



(3.2) 



V c (k — k') = -8ttk 



Afi 2 



[—(k - k') 2 + /i 2 ] 2 [—(k - k') 2 + n 



21:! 



(3.3) 



V c (r, T) = nr exp[-/x(T)r] , (3.4) 



" (T) = 1 - 0.7?T/T C ) 2 ' (3 ' 5) 



, s « 1 - 0.7(T/T C ) 2 1 , , 

= ^ W CJ J • (3.6) 

/i e 
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Fig. 5 A comparison of quenched (open symbols) and unquenched (filled symbols) results 
for the interquark potential at finite temperature. The dotted line is the zero temperature 
quenched potential. Here, the symbols for T = 0.80T C [open triangle], T = 0.88T C [open 
circle], T = 0.94T C [open square], represent the quenched results. The results with dynamical 
fermions are given at T = 0.68T C [solid downward-pointing triangle], T = 0.80T C [solid 
upward-pointing triangle], T = 0.88T C [solid circle], and T = 0.94T C [solid square]. 
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Fig. 6 The potential V c (r, T) is shown for T/T c = [solid line], T/T c = 0.4 [dotted 
line], T/T c = 0.6 [dashed line], T/T c = 0.8 [dash-dot line], T/T c = 0.9 [short dashes], T/T c = 
1.0 [dash-dot-dot line]. Here, V c (r,T) = «rexp[-^(T)r], with /i(T) = 0.01GeV/[l - 
0.7(T/T C ) 2 ]. 
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Fig. 7 Temperature dependent constituent mass values, m u (T) and m s (T), calculated 
using the equation below, are shown. Here mP u = 0.0055 GeV, m? s = 0.120 GeV, and 
G(T) = 5.691 [1 - 0.17(T/T C )], if we use Klevansky's notation. 
m(T) = m° + 2Gs{T)N c ^ J Q A dpf- tanh(±/3£ p ). 
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Fig. 8 Values of m u (T) are shown. The dashed curve is calculated with m° = 5.50 MeV. 
Here, G{T) =G[1- 0.17 (T/T c )\, with G = 5.691 GeV" 2 and T c = 0.150 GeV. The solid 
curve is calculated with the same value of G(T) and T c , but with m° = 0. From the solid 
curve, we see that chiral symmetry is restored at T = 0.136 GeV when m° = 0. 
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Fig. 9 The mass values of the pionic states calculated in this work with G n (T) = 
13.49[1 - 0.17T/TJ GeV, G V {T) = 11.46[1 - 0.17T/TJ GeV. The value of the pion mass is 
0.223 GeV at T/T c = 0.90, where m u (T) = 0.102 GeV and m s (T) = 0.449 GeV. The pion 
is bound up to T/T c = 0.94, but is absent beyond that value. 
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IV. CALCULATION OF HADRONIC CURRENT CORRELATION FUNCTIONS 
AT FINITE TEMPERATURE 

Bing He, Hu Li, C. M. Shakin, and Qing Sun, Physical Review C 67, 065203 (2003) 



Fig. 10 The upper figure represents the basic polarization diagram of the NJL model 
in which the lines represent a constituent quark and a constituent antiquark. The lower 
figure shows a confinement vertex [filled triangular region] used in our earlier work. For the 
present work we neglect confinement for T > 1.2 T c , with T c = 150 MeV. 

For ease of reference, we present a discussion of our calculation of hadronic current cor- 
relators. The procedure we adopt is based upon the real-time finite-temperature formalism, 
in which the imaginary part of the polarization function may be calculated. Then, the real 
part of the function is obtained using a dispersion relation. The result we need for this 
work has been already given in the work of Kobes and Semenoff. (The quark momentum is 
k and the antiquark momentum is k — P. We will adopt that notation in this section for 
ease of reference to the results presented in the work of Kobes and Semenoff.) We write the 
imaginary part of the scalar polarization function as 








{(1 - ni (k) - n 2 (k))5(P° - E,(k) - E 2 {k)) 



(4.1) 



-MA;) - n 2 (k))5(P° + E 1 (k) - E 2 {k)) 
-M*0 - m(k))S(P° - E ± (k) + E 2 {k)) 
(1 - m(k) - n 2 (k))5(P° + E 1 (k) + E 2 (k))} . 
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Here, E x {k) = [ k 2 + m\{T) ] 1//2 . We have included a Gaussian regulator, exp[— k 2 /a 2 ], with 
a = 0.605 GeV, which is the same as that used in most of our applications of the NJL model 
in the calculation of meson properties. We also note that 

^) = e ^l) + ! . ( 4 - 2 ) 

and 

n2 ^ = e ^(l) + x • ( 4 - 3 ) 

For the calculation of the imaginary part of the polarization function, we may put k 2 = 
m 2 (T) and (k — P) 2 = m\(T), since in that calculation the quark and antiquark are on- 
mass-shell. In Eq. (14. lj) the factor /3 s arises from a trace involving Dirac matrices, such 
that 

(3 S = -Ti[^ + mi )^-f + m 2 )] (4.4) 
= 2P 2 -2(mi +m 2 ) 2 , (4.5) 

where mi and m 2 depend upon temperature. In the frame where P = 0, and in the case 
mi — m 2i we have (3g = 2P 2 (1 ~ <im 2 / Pq). For the scalar case, with m\ = m 2 , we find 

W S (P*,T) = ^ (l - ^)' V2 e- J2/ " 2 [l - am(*)], (4-6) 



where 



k 2 = tt-m 2 (T). (4.7) 

For pseudoscalar mesons, we replace by 

/3 P = -Trf^^ + mO^^-f + m 2 )] (4.8) 
= 2P 2 - 2(mi - m 2 ) 2 , (4.9) 

which for rri\ = m 2 is f3p = 2P 2 in the frame where P = 0. We find, for the 7r mesons, 

hnMf.T) = ^ (l - ^) V ^- SI/ " 3 [l - *».(*>], (4-10) 

where A; 2 = Pq/4 — m 2 (T), as above. Thus, we see that, relative to the scalar case, the 
phase space factor has an exponent of 1/2 corresponding to a s-wave amplitude. For the 
scalars, the exponent of the phase-space factor is 3/2. 
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For a study of vector mesons we consider 

= Tr[ 7 ^ + mi ) lv {$ - f + m 2 ) ] , (4.11) 

and calculate 

= 4[P 2 - m?-m^ + 4mim 2 ], (4.12) 

which, in the equal- mass case, is equal to 4P 2 + 8m 2 (T), when mi = m 2 and P = 0. This 
result will be needed when we calculate the correlator of vector currents in the next section. 
Note that, for the elevated temperatures considered in this work, m u (T) = md{T) is quite 
small, so that 4Pq + 8m 2 (T) can be approximated by 4P 2 , when we consider the vector 
current correlation functions. In that case, we have 

Im J V (P 2 ,T) ~ |lm J P (P 2 ,T). (4.13) 
At this point it is useful to define functions that do not contain that Gaussian regulator: 

Im J P (P 2 , T) = ( 1 _ V2 [ 1 - 2ni(k) } , (4.14) 



o 



and 

Im J V {P\ T) = H^S _ V2 [! - 2n 1 (Jb) ] , (4.15) 

We need to use a twice-subtracted dispersion relation to obtain ReJp(P 2 ,T), or 
ReJ y (P 2 ,T). For example, 
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Re J P (P 2 , T) = Re J P (0, T) + Re J P (P 2 , T) - Re J P (0, T) ] + (4.16) 

o 

P 2 (P 2 -P 2 ) / A2 ^_ Im J P ( S ,T) 



74 m 2 (r) S(P 2 - S)(P 2 - S) ' 

where A 2 can be quite large, since the integral over the imaginary part of the polarization 
function is now convergent. We may introduce Jp(P 2 ,T) and Jy(P 2 ,T) as complex func- 
tions, since we now have both the real and imaginary parts of these functions. We note that 
the construction of either Re Jp(P 2 ,T), or Re Jy(P 2 ,T), by means of a dispersion relation 
does not require a subtraction. We use these functions to define the complex functions 
J P {P 2 ,T) and J V {P 2 ,T). 
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In order to make use of Eq. (j4.16|) . we need to specify Jp(0) and J P (P 2 ). We found it 
useful to take P 2 = -1.0 GeV 2 and to put J P (0) = J P (0) and J P (P 2 ) = J p (Pq). The 
quantities Jy(0) and Jy(P 2 ) are determined in an analogous function. This procedure in 
which we fix the behavior of a function such as ReJy(P 2 ) or ReJy(P 2 ) is quite analogous to 
the procedure used in our earlier work. In that work we made use of dispersion relations to 
construct a continuous vector-isovector current correlation function which had the correct 
perturbative behavior for large P 2 — > — oo and also described that low-energy resonance 
present in the correlator due to the excitation of the p meson. In our earlier work the 
NJL model was shown to provide a quite satisfactory description of the low-energy resonant 
behavior of the vector-isovector correlation function. 

We now consider the calculation of temperature-dependent hadronic current correlation 
functions. The general form of the correlator is a transform of a time-ordered product of 
currents, 

iC(P 2 ,T) = J d 4 xe lP - x « T(j(x)j{0)) » , (4.17) 

where the double bracket is a reminder that we are considering the finite temperature case. 

For the study of pseudoscalar states, we may consider currents of the form jp t i{x) = 
q(x)i'j5\ l q(x), where, in the case of the 7r mesons, i — 1,2 and 3. For the study of scalar- 
isoscalar mesons, we introduce js,i{x) — q{x)X l q(x), where i = for the flavor-singlet current 
and i = 8 for the flavor-octet current. 

In the case of the pseudoscalar-isovector mesons, the correlator may be expressed in terms 
of the basic vacuum polarization function of the NJL model, J P (P 2 ,T). Thus, 

where Gp(T) is the coupling constant appropriate for our study of 7r mesons. We have 
found Gp(T) = 13.49 GeV~ 2 by fitting the pion mass in a calculation made at T = 0, with 
m u = irid = 0.364 GeV. The result given in Eq. 1)4.18)1 is only expected to be useful for small 
P 2 , since the Gaussian regulator strongly modifies the large P 2 behavior. Therefore, we 
suggest that the following form is useful, if we are to consider the larger values of P 2 . 



C P {P\T) 
P 2 



Jp{P\T) 
P 2 



' (4.19) 



l-Gp(T)Jp(P 2 ,T)- 
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(As usual, we put P = 0.) This form has two important features. At large Pq, 
ImCp(Po, T)/Pq is a constant, since Im Jp(P 2 ,T) is proportional to Pq. Further, the 
denominator of Eq. ()4.19|) goes to 1 for large Pq. On the other hand, at small P 2 , the 
denominator is capable of describing resonant enhancement of the correlation function. As 
we will see, the results obtained when Eq. (14.19)1 is used appear quite satisfactory. 

For a study of the vector-isovector correlators, we introduce conserved vector currents 
jnA x ) = Q( x )l^i ( l( x ) with i=l, 2 and 3. In this case we define 

pfj.pi 



and 



0?{P\T)=(g 



p2 



pfipv 



Jv{P\T) 



(4.20) 



p2 



C V {P\T) 



(4.21) 



taking into account the fact that the current j^i(x) is conserved. (Note that Eqs. (|4.20|) and 
(|4.21|) are valid for zero temperature. However, we still use that form at finite temperature 
for convenience.) We may then use the fact that 



J v {P\T)= l -g, u J^{P\T) 



(4.22) 



and 



rp 2 + 2m 2 (T) 



Im J V (P 2 ,T) 



4tt 



-ImJ P (P 2 ,T) 



4m 2 (T) 



p2 




1/2 



e - k2 / a2 [l-2n 1 {k)} (4.23) 

(4.24) 



(See Eq. (J4.7)) for the specification of k = \k\.) We then have 

1 



C V {P\T) = J V {P\T) 



l-G v (T)J v (P\T) 



where we have introduced 



ImJ y (P 2 ,T) 



P 2 + 2m 2 (T) 

47T 



-ImJp(P 2 ,T). 



4m 2 (T)\ 1/2 



P 2 



2 ni (k) 



In the literature, u is used instead of Pq. We may define the spectral functions 

ay(uj, T) = — Im CV(cj, T) , 

7T 



(4.25) 



(4.26) 
(4.27) 



(4.28) 
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and 

a P {u,T) = -lmC P {u,T) , (4.29) 

7T 

Since different conventions are used in the literature, we may use the notation Wp(uj,T) 
and <5V(u;,T) for the spectral functions given there. We have the following relations: 

a P (uj,T) = ap(uj,T), (4.30) 

and 

a v (uj,T) 3 

2 = l a v{u, T) , (4.31) 

where the factor 3/4 arises because there is a division by 4 in the literatures, while we have 
divided by 3. 

V. CALCULATION OF THE MOMENTUM DEPENDENCE OF HADRONIC 
CURRENT CORRELATION FUNCTIONS AT FINITE TEMPERATURE 



Xiangdong Li, Hu Li, C. M. Shakin, Qing Sun and Huangsheng Wang, nucl-th/0405081 



We have calculated spectral functions associated with hadronic current correlation func- 
tions for vector currents at finite temperature. We made use of a model with chiral sym- 
metry, temperature-dependent coupling constants and temperature-dependent momentum 
cutoff parameters. Our model has two parameters which are used to fix the magnitude 
and position of the large peak seen in the spectral functions. In our earlier work, good fits 
were obtained for the spectral functions that were extracted from lattice data by means 
of the maximum entropy method (MEM). In the present work we extend our calculations 
and provide values for the three-momentum dependence of the vector correlation function 
at T = 1.5 T c . These results are used to obtain the correlation function in coordinate 
space, which is usually parametrized in terms of a screening mass. Our results for the 
three-momentum dependence of the spectral functions are similar to those found in a recent 
lattice QCD calculation for charmonium [S. Datta, F. Karsch, P. Petreczky and I. Wet- 



zorke, hep-lat/0312037|. For a limited range we find the exponential behavior in coordinate 
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space that is usually obtained for the spectral function for T > T c and which allows for the 
definition of a screening mass. 




Fig. 11 The spectral functions a/uo 2 for pseudoscalar states obtained by MEM are 
shown. The solid line is for T/T c = 1.5 and the dashed line is for T/T c = 3.0. The second 
peak is a lattice artifact. 



0.18-, 




co/T 

Fig. 12 The spectral functions a/uo 2 for vector states obtained by MEM are shown. The 
second peak is a lattice artifact. 
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Fig. 13 The imaginary part of the correlator a(u)/ou 2 is shown for various values of \P\ 
as a function of uj 1 . Starting with the topmost curve the values of \P\ in GeV units are 0.10, 
0.30, 0.50, 0.70, 0.90, 1.10, 1.30, 1.50, 1.70, 1.90 and 2.10. Here we have used G s = 1.2 
GeV" 2 and k max = 1.22 GeV. 
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Fig. 14 The correlation function C{z) is shown. The dotted line represents a fit using 
an exponential function. 

Here, C{z) = \ dP z e %p * z J °° du a{uJ ^ ' Pz) . We may also use the form C{z) = 
\ /_~ dP z e*" J °° dP* ^hhl . 
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VI. QUARK PROPAGATION IN THE QUARK-GLUON PLASMA 

Xiangdong Li, Hu Li, C. M. Shakin and Qing Sun, Physical Review C, 69, 065201 

(2004) 

It has recently been suggested that the quark-gluon plasma formed in heavy-ion collisions 
behaves as a nearly ideal fluid. That behavior may be understood if the quark and antiquark 
mean-free- paths are very small in the system, leading to a "sticky molasses" description 
of the plasma, as advocated by the Stony Brook group. This behavior may be traced to 
the fact that there are relatively low-energy qq resonance states in the plasma leading to 
very large scattering lengths for the quarks. These resonances have been found in lattice 
simulation of QCD using the maximum entropy method (MEM). We have used a chiral 
quark model, which provides a simple representation of effects due to instanton dynamics, 
to study the resonances obtained using the MEM scheme. In the present work we use our 
model to study the optical potential of a quark in the quark-gluon plasma and calculate 
the quark mean-free-path. Our results represent a specific example of the dynamics of the 
plasma as described by the Stony Brook group. 




P 2 (GeV 2 ) 



Fig. 15 Values of t(P 2 ,p2) are shown for various values of the quark momentum \p 2 \. 
Starting with the uppermost curve, the \p 2 \ values in GeV units are 0.01, 0.03, 0.05, 0.07, 
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0.09, 0.11, 0.13, 0.15, 0.17, 0.19, 0.21, 0.23, 0.25, 0.27, 0.29 and 0.31. (For large P 2 , we have 
t(P 2 ,p 2 ) ^ (l/nP 2 )G.) Here P 2 = (pi + P2) 2 , where p\ is the antiquark momentum. 




p, (GeV) 



Fig. 16 Values of n(pi) are shown for /i — 1.1 GeV (dotted curve), \x = 1.3 GeV (dashed 
curve) and /1 = 1.5 GeV (solid curve). Here T = 1.5 T c with T c = 270 MeV. 
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Fig. 17 The imaginary part of the quark optical potential is shown for fi = 1.1 GeV 
(dotted curve), fx — 1.3 GeV (dashed curve) and fx — 1.5 GeV (solid curve). (We recall that 
the nucleon-nucleus imaginary optical potential is about 0.01 GeV in magnitude.) 
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VII. CALCULATION OF SCREENING MASSES IN A CHIRAL QUARK MODEL 

Xiangdong Li, Hu Li, C. M. Shakin and Qing Sun, nucl-th/0405035 



We consider a simple model for the coordinate-space vacuum polarization function which 
is often parametrized in terms of a screening mass. We discuss the circumstances in which 
the value m sc = ttT is obtained for the screening mass. In the model considered here, that 
result is obtained when the momenta in the relevant vacuum polarization integral are small 
with respect to the first Matsubara frequency. 

In order to present our results in the simplest form, we consider only the scalar interaction 
proportional to (qq) 2 . We also extend the definition of a(u, T) to include a dependence upon 
the total moment of the quark and antiquark appearing in the polarization integral. Thus 
we consider the imaginary part of the correlator, a(u, P). Since we place P along the z-axis 
this quantity may be written as cr(u, 0, 0, P z ). In this work we will present our result for the 
coordinate-dependent correlator C(z) which is proportional to the correlator, 

2 J-oo Jo ^ 

We may also use the form 

C(z) = \ |°° dP z e iP * z J™ dP 2 . (7.2) 

We have made a study of the screening mass in a simple model in order to understand 
the origin of exponential behavior for the correlator. We consider the Matsubara formalism 
and note that the quark propagator may be written, with f3 = 1/T, 

Sp(k,uj n ) = => . (7.3) 

(2n + l)% 2 //3 2 + k 2 + M 2 

For bosons the vacuum polarization function is given as, 

it/ — > o\ g 2 ^d?k 1 1 



P 2 

We modify the last equation to refer to fermions. In this case the Matsubara frequencies 

are 

(2n + l)n 

UJ n = (7.5) 
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and we have 



<r ... f d 3 k 





7 V + tt//3) + t" • ("fc + ~p )) 















, (7.6) 



if we keep only the first term in the sum, where ujq — ir/(3. As a next step we drop p°, so 
that we have 



n(?,o) = |> ' (H ' 



(VV,3 + ~i ■ ~k ) (tV/? + 7" • (k + -p )) 



2/? 7 (2?r) 



7T 



7T 



73) - (/- +" 



(7.7) 



We then take p along the z axis and write II (p z ) = II( p , 0). We define 

C(z) = y rfp 2 e^ 2 n(p 2 ). (7.8) 

In our calculation we replace g 2 /2[3 by unity and use a sharp cutoff so that | k \ < k max . 
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Fig. 18 The function C(z) is shown for a sharp cutoff of k max = 0.1 GeV. The dotted 
line represents an exponential fit to the curve using m sc = 1.23 GeV. (We recall that irT is 
equal to 1.27 GeV.) 
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Fig. 19 The function C(z) of is shown for a sharp cutoff of k ma x = 0.4 GeV. The dotted 
line represents an exponential fit to the curve using m sc = 0.961 GeV. (We recall that nT is 
equal to 1.27 GeV.) 
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